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ABSTRACT 

We present results from a comprehensive imaging survey of 70 radio galaxies at redshifts 1 < 
z < 5.2 using all three cameras onboard the Spitzer Space Telescope. The resulting spectral energy 
distributions unambiguously show a stellar population in 46 sources and hot dust emission associated 
with the active nucleus in 59. Using a new restframe S3 ^m/ Si.6/xm versus S5 urn/ S3 /im criterion, we 
identify 42 sources where the restframe 1.6 /xm emission from the stellar population can be measured. 
For these radio galaxies, the median stellar mass is high, 2 x 10 u Mq, and remarkably constant 
within the range 1 < z < 3. At z > 3, there is tentative evidence for a factor of two decrease in stellar 
mass. This suggests that radio galaxies have assembled the bulk of their stellar mass by z ~ 3, but 
confirmation by more detailed decomposition of stellar and AGN emission is needed. 

The restframe 500 MHz radio luminosities are only marginally correlated with stellar mass but are 
strongly correlated with the restframe 5 /im hot dust luminosity. This suggests that the radio galaxies 
have a large range of Eddington ratios. We also present new Very Large Array 4.86 and 8.46 GHz 
imaging of 14 radio galaxies and find that radio core dominance — an indicator of jet orientation — 
is strongly correlated with hot dust luminosity. While all of our targets were selected as narrow-lined, 
type 2 AGNs, this result can be understood in the context of orientation-dependent models if there 
is a continuous distribution of orientations from obscured type 2 to unobscured type 1 AGNs rather 
than a clear dichotomy. Finally, four radio galaxies have nearby (<6") companions whose mid-IR 
colors are suggestive of their being AGNs. This may indicate an association between radio galaxy 
activity and major mergers. 

Subject headings: galaxies: active — galaxies: evolution — galaxies: high-redshift — radio continuum: 
galaxies 



1. INTRODUCTION 

Across cosmic time, powerful radio sources are 
robust beacons of the most massive galaxies in the 
Universe. At low redshift, this has been known since 
the first visible counterparts of cxtragalactic radio 
sources were shown to be gian t elliptical (gE and cD) 
galaxies (Matthew s et al.l [1964) . In the more distant 
Universe, indirect evidence for this correlation initially 
came from a variety of observations, including the 
near-IR r 1 / 4 light profiles of high redshift radio galaxies 
(HzRGs) at 1 < z < 2 in Hubb l e Space Te l escope 
images (e.g., iPentericci etafl 120001: iZirm et al.l [2003) 
and the tendency for HzRGs to resid e in moderately 
rich Jproto-) cluster environ ments (e.g., iVenemans et al.l 
Galametz et al.l [20 10). The most direct evidence 
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comes from the remarkably tight scatter of the observed 
if- band magn itudes of HzRGs in the Hubble K-z 
diagram (e.g.. iLillv fc Longairl Il984t lEales et al.l 119971: 
Best et all 119981: Ivan Breugel et al.l 119981: Uarvis et al.l 
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Bryant et al.l 12009). Assuming that the observed- frame 



Willott e t alJ 200; 
iBrookes et atl [200! 



if-band light is dominated by emission from old stellar 
populations, these studies put HzRGs at the top end 
of the stellar mass function out to z — 5.2. However, 
this assumption is uncertain due to large band-shifting 
effects and remaining contributions from the AGN and 
young stellar populations, which undoubtedly contribute 
to the observed if -band magnitudes particularly at the 
highest redshifts where the observed if-band measures 
rest-frame UV light. As Seymour et al. (2007, hereafter 
S07) argued in the precursor to this paper, the only 
solution to minimize these uncertainties is to observe 
at longer wavelengths with the Spitzer Space Telescope. 
Such observations avoid k-correction effects by consis- 
tently observing the same restframe wavelengths where 
the old stellar population peaks. S07 argued that the 
restframe iZ-band luminosity (corrected for non-stellar 
emission) is the most efficacious band for deriving stellar 
masses because it corresp onds to the minimum in the 
opacity of the H~ ion fe.g.. lJohn"lll988H and is associated 
with a bump in the stellar population SEP for almost all 

Simpson fc Eisenhardl U99l 
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stellar populations 
see S07 for details) 

The initial observations for our program, the Spitzer 
High- Redshift Radio Galaxy or SHzRG project, were ob- 
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tained during Spitzer Cycle 1 (Program ID 3329) and 
involved mid-IR (A > 3.5 /im) imaging of 70 HzRGs at 
1 < z < 5.2. S07 reported the sample selection, data 
processing, and initial results from the full survey. Us- 
ing imaging from all three instruments on Spitzer, S07 
decomposed the restframe visible to infrared spectral en- 
ergy distributions (SEDs) of HzRGs into stellar, AGN, 
and dust components and determined the contribution 
of host galaxy stellar emission at restframe i7-band to 
derive stellar masses. S07 found that >60% of restframe 
ff-band light is from stars for the majority of well-studied 
HzRGs. As expected from unified models of AGNs, the 
fraction of restframe ii-band luminosity due to stars is 
not correlated with redshift, radio luminosity, or rest- 
frame mid-IR (5 /im) luminosity. In addition, while the 
stellar 77-band luminosity was not found to vary with 
stellar fraction, the total 77-band luminosity was found 
to anti-correlate with stellar fraction, as expected if the 
underlying HzRG hosts comprise a homogeneous popu- 
lation. From a comparison with predicted luminosities 
from stellar evolution models, S07 found that the HzRG 
restframe H-band stellar luminosities imply host stellar 
masses of 10 11 -10 115 M for HzRGs, even at the high- 
est redshifts, implying an early formation epoch for these 
massive galaxies. However, with a large fraction of the 
sources poorly sampled longward of observed 8 /zm, many 
of the results of S07 were based on upper limits or on the 
subset of galaxies with well-sampled SEDs. 

In addition to the comprehensive results presented 
by S07, our SHzRG observations have been used to 
comp lement samples of lower redshift radio galaxies 
(e.g., lHaas et al.l 12008ft and to explore HzRG environ- 



ments (e.g.. iGalametz et al.l 120091: i Dohertv et al.l 120101: 
Galametz et al.l 120101 : iKuiper et al.l 12010b iTanaka et al.l 
2010ft . The legacy value of the SHzRG program is 
also seen in the numerous studies of individual HzRGs 
which have drawn upon the SHzRG data set — e.g., 
Stern et al. (2006; LBDS 53W091), Villar-Martin et 
al. (2006; MRC 2104-242), Broderick et al. (2007; 
PKS 0529-549), Greve et al. (2007; 4C 41.17), Ivi- 
son et al. (2008; 4C 60.07), Nesvadba et al. (2009; 
USS 0828+193), Smith et al. (2010; TXS J1908+7220), 
and De Breuck et al. (in prep.; 4C 23.56). Follow-up 
mid-IR spectroscopy of two SHz RG targets (4C 23.5 6 
and 4C 72.26) is also reported bv lSevmour et al.l (|2008ft . 

This paper reports the results of a Cycle 4 GTO pro- 
gram (Program ID 40093) to complete the Spitzer ob- 
servations of the SHzRG sample. The observations pre- 
sented here entail complete six-band 3.6 to 24 /zm mid- 
IR imaging photometry for the entire SHzRG sample 
using all three instruments on-board Spitzer. In addi- 
tion, new radio observations complete the radio mor- 
phological parameters for our sample. Throughout we 
assume a concordance model of Universe expansion with 
ti M = 1-^A = 0.3, tt = 1, andiJo = 70 kms" 1 Mpc -1 . 
Infrared luminosities presented in this paper are of the 
form vLv/Lq, where Lq = 3.9 x 10 26 W. Following com- 
mon practice, radio luminosities at 500 MHz and 3 GHz 
are expressed as monochromatic luminosity densities. 



2. DATA 
2.1. Sample 



The SHzRG sample was described in detail by S07 and 
is listed in Table [1] In short, we selected 69 powerful 
radio galaxieifl from the literature covering a uniform 
range in both redshift (requiring z > 1) and 3 GHz rest- 
frame luminosity. In practice, all z > 3.5 radio galax- 
ies known in 2004 February were included in our sam- 
ple. As no low-power radio sources are known as yet 
at the highest redshifts, we unfortunately are unable 
to study radio power dependence at z > 3. At lower 
redshifts, we preferentially selected those sources with 
the most supporting data. During Cycle 1, we inadver- 
tently observed MRC 0211-256 (z = 1.300) instead of 
TXS 0211-122 (z = 2.340). As the latter is a well- 
studied radi o galaxy with suppor t ing Hubble Space Tele- 
scope data (Ivan Oiik et all 119971 : iPentericci et al.l 120011: 
iVernet et al.ll200lD . we have added this source to our 
sample in Cycle 4, which thus now contains 70 radio 
galaxies at 1.0 < z < 5.2. 
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Figure 1. Restframe 3 GHz luminosity plotted against restframe 
500 MHz luminosity. The low scatter shows that the contribution 
by Doppler-boosted radio emission in our sample is nearly always 
negligible. The solid line shows spectral index a = —1.0, typical of 
the radio galaxies in our sample. 

We initially calculated the radio luminosities at a rest- 
frame frequency of 3 GHz in order to be able to inter- 
polate between the all-sky, low-frequency (327/352, 365 
and 1400 MHz) ra dio surveys available at the time of our 
Cycle 1 proposal (iRengelink et al.lll997t iDouglas et al.l 
[19961: ICondon et al .||1998l). Since lhen, the 74 MHz VLA 
Low Frequency Survey (VLSS; iCohen et al.l l200"7ft has 
been com pleted, covering the entire sky at 5 > —30°. 
Following iMilevfe De Breuckl (|2008ft . we therefore now 
also calculate restframe 500 MHz luminosities for our 
sample using the VLSS and the NRAO VLA Sky Survey 
(NVSS; ICondon et all [HUH- Table [2] lists both the 

9 We define radio galaxies as galaxies (i.e., not type 1 unob- 
scured quasars) with restframe 3 GHz luminosities greater than 
10 26 WHz" 1 . 

10 For the few sources not detected in the VLSS, we 
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500 MHz and 3 GHz restframe radio luminosities. The 
median flux ratio between the two frequencies is 6.3 for 
our sample (see Fig. [1]), corresponding to a spectral in- 
dex a = —1.0 (S v oc v a ). This steeper-than-average 
spectral index reflects the steep spectrum selection cri- 
teria adopted in most of the parent samples from which 
our targets were drawn. 

Figure [2] revisits the radio power-redshift plane used 
in the definition of our sample. Compared to the 3 GHz 
luminosity versus redshift plot (bottom panel), there is 
a slightly larger degeneracy with Lsoomhz (top panel), 
particularly at z > 3.5 and z < 2. Our sample still 
contains sources with log(LsooMHz) ~ 28.5 throughout 
the redshift range 1 < z < 4.5 but covers only ~1 dex 
at a given redshift. As argued by S07, we do not expect 
Doppler-boosted emission to significantly contribute in 
our sample, even at restframe 3 GHz. However, in this 
paper we will use the 500 MHz rather than the 3 GHz 
luminosities because emission at the low er frequency is 
expec ted to be more nearly isotropic (e.g.. lBlundell et all 
1998). 
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Figure 2. Restframe 500 MHz (top) and 3 GHz (bottom) lumi- 
nosity plotted against redshift. Redshift and radio power are more 
degenerate at 500 MHz. 



2.2. Spitzer observations 

Our Cycle 1 SHzRG observations used a combination 
of the imaging modes for all three instruments onboard 
Spitzer. We observed (i) all 69 sources from our (initial) 
sample in all four band^3 — 3.6, 4. 5, 5.8 and 8.0 fim — of 
the Infrared Array Camera (IRAC: lFazio et al.ll2004[ ). (ii) 
all 46 sources at z > 2 using the 16 /im pea k-up imaging 
came ra of the InfraRed Spectrograph (IRS; iHouck et al.l 
l200l . and (iii) all 26 sources with "low" predicted mid- 
IR background (SWum < 20 MJy sr _1 ) in all three bands 
— 24, 70, and 160 — of the Multiband Im aging Pho- 
tometer for Spitzer (MIPS; IRieke et al.ll200l . The Cy- 
cle 4 program completes the full six-band 3.6 to 24 /im 

used the 325 MHz W esterbork Northern Sky Sur vey (WENSS; 
Rcngclink et al. 1997) or the 365 MHz Texas survey (Douglas et al. 



1996). 



With the exception of 3C 65, which was only observed in the 
3.6 and 5.8 fim bands (see S07). As the source is dominated by hot 



dust emission at A b s 
8.0 fim during Cycle 4 



5.8 ^m, we did not re-observe it at 4.5 and 



photometry. Table Q] summarizes all Spitzer imaging ob- 
servations of our sample, including those previously pub- 
lished by S07. A handful of sources have significantly 
deeper Spitzer exposures in one or more bands (e.g., deep 
IRAC observations of 4C 41.17; IRAC and MIPS obser- 
vations of B2 0902+34). These were obtained by other 
programs (GO and GTO) and are included here as part 
of the SHzRG sample. We now describe the Cycle 4 
SHzRG program in more detail. Table [3] lists the full 
six-band photometry for our sample, and Appendix A 
shows the individual SEDs. 

IRAC: We observed only TXS 0211-122 (see §2.1) 
with IRAC during Cycle 4. The observations were iden- 
tical to our Cycle 1 program, consisting of four dithered 
exposures in each of the channels. Data reduction fol- 
lowed the methods described by S07 except that we used 
IRAC pipeline version 16.1 instead of 13.2. As there are 
no calibration differences between the pipelines, we did 
not re-reduce the other IRAC data. 

IRS: We observed the 24 z < 2 objects using the 16 fim 
peak-up imaging camera of IRS. During Cycle 1, this ob- 
serving mode was not officially supported, and we used 
an experimental double nodding sequence. For the new 
Cycle 4 observations, we switched to the more efficient 
peak-up astronomical observing template consisting of 
five random positions of 30 s ramp duration each. To en- 
sure a fair comparison between our Cycle 1 and Cycle 4 
photometry, we re-obtained all our data from the Spitzer 
archive (pipeline version S16.0.0) and re-mosaiced the 
data with M0PEX. Since the Cycle 1 data consisted of just 
two exposures, we mosaiced the difference images after 
subtracting a median to remove the background (includ- 
ing residual images in a few cases). The final images 
were resampled on a grid with l'/8xl'.'8 pixels. We ex- 
tracted the fluxes in circular apertures with a 6-pixel 
radius and a background subtraction annulus between 6 
and 10.7 pixels. We applied an aperture correction of 
1.154 to our photometry to obtain total fluxeJ^l. The 
new fluxes listed in Table [3] are on average 50 % higher 
than the one published by S07 due to the improved cali- 
bration. 

MIPS: Our Cycle 1 observations gave a low detection 
rate in the MIPS 70 fim imaging (five out of 24: 21%) 
and no detections at 160 fim. The 44 HzRGs not ob- 
served with MIPS during Cycle 1 were omitted because 
of the larger predicted background emission (S^^m > 
20MJysr _1 , see S07) at their positions. Because of the 
low detection rate at 70 and 160 fim and the higher ex- 
pected background emission due to zodiacal and Galactic 
cirrus emission, an observing request aimed at detect- 
ing our complete sample at these wavelength was con- 
sidered impractical. We therefore decided to concentrate 
on 24 /im imaging during Cycle 4. To compensate for 
the higher expected sky background, the new 24 fim ob- 
servations consisted of two cycles of 30 s exposures each 
(totaling 28 individual images compared to a single cy- 
cle, i.e., 14 images, in our Cycle 1 observations). One 
radio galaxy, MRC 0156—252, was not observed before 
the end of the Spitzer cryogenic mission. To ensure a fair 
comparison of flux densities between Cycles 1 and 4, we 
re-obtained the data from the archive (pipeline version 
SI 6. 1.0) and mosaiced and derived photometry in the 

12 See http: / /ssc. spitzer. caltech.edu/irs/irsinstrumenthandbook/45/ 



4 



De Breuck et al. 




Figure 3. Restframe luminosity density of the Spit zer data in our sam ple. The red triangles show the average 3CRR radio galaxy SED 
and the blue diamonds the average quasar SED from Lcipski et al. (2010), both based on objects with 1 < z < 1.4 a nd normalized to radio 
flux density. The green solid line is the mean SED of all SDSS (visible color selected, radio quiet) type-1 AGNs from Richards ct al. (2006) 
normalized to 0.08 mjy at 2500 A. 



same fashion as S07. Sources with predicted sky back- 
ground levels S24 fim > 40 MJy,sr _1 are affected by strong 
remaining gradients in the combined images. While such 
fields are not appropriate to study the environment of 
the HzRGs (Mayo et al., in prep.), the images can still 
be used to obtain reliable photometry for the central ra- 
dio galaxies. The one exception is 5C 7.269, the source 
with the highest observed background level in our sam- 
ple (5 , 24 M m = 65MJy,sr _1 ), which saturated the MIPS 
detector. 

2.3. New radio data 

Of the 70 sources in our SHzRG sample, 49 already 
had high resolution (0'.'5) radi o maps available f rom 
the Very Large Array (VLA; INapier et al.l 1 19831) or 
the Australia Telescope C ompact Array fATCA; iFraterl 
iBrooks fc Whiteoak 1992f ). Most of these observations 
are co ntaine d in the compend i a publ ished bv lCarilli et al.l 
(|1997|) and iPentericci etaTl 1)20000 . To complete the 
radio observations for our sample, we observed the re- 
maining 21 sources with the VLA in the A configura- 
tion between 2006 February 20 and 2006 March 16 (ob- 
serving project ID AD520). F or consistency, we used 
the same observationa l setup as iCarilli et al.l (|1997| ) and 
IPentericci et al.l (|2000l ); i.e., 10 to 40min of snapshot ob- 
servations in C-band (4.86 GHz) and 20 to 80min in X- 
band (8.46 GHz). We followed the standard calibration 



and imaging reduction steps in the Astronomical Imag- 
ing Processing System (AXPS) and the Common As- 
tronomy Software Applications (CASA). 

Table [5] summarizes the archival and new radio data for 
our sample. We classified the radio structures into three 
classes: S — single component (i.e., spatially unresolved 
sources); D = double component sources without a core 
detection; and T = core-detected, resolved sources, which 
are generally triple component sources with a core and 
two radio lobe s. For all sources, w e measured the largest 
angular size 9 (jCarilli et al.lll997l) of the radio structure. 
For sources with a core detection, we list the core fraction 
at a restframe frequency of 20 GHz, CF20 (following the 
definition of Carilli et al. 1997). Appendix B shows 
contour plots of the new radio data. 

3. MID-IR SPECTRAL ENERGY DISTRIBUTIONS 

3.1. Mean HzRG SED 

As our Spitzer survey presents the largest collection 
of mid-IR photometry of targeted, high-redshift type 2 
AGNs, we first examine the general trends in their SEDs. 
Figure |3] shows the restframe infrared luminosity den- 
sities from our Spitzer photometry. There is a scatter 
of ^2 dex at any given rest wavelength with a notable 
exception around A rcs t = 1.6 //m, where the scatter de- 
creases to 1.5 dex. This wavelength corresponds to the 
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peak in old stellar populations (e.g., Sawicki 2002), and 
the decreased scatter suggests that the stellar population 
dominates the SEDs at this wavelength. The 1.6 /im peak 
is most pronounced in the lowest luminosity objects and 
tends to become less obvious in sources with greater mid- 
IR luminosity. Almost all the SEDs show a sharp rise at 
A rest > 2 /im due to the onset of hot dust emission at 
longer wavelengths. 

The avera g e rad io galaxy SED as derived by 
iLeipski et all pOlCh is a good match to our sample, 
though the most luminous sources have SEDs compat- 
ible with either of the two quasar SEDs shown in Fig- 
ure [31 The radio galaxy points have a steeper rise be- 
tween restframe ^2 and ~8 /im than the quasar com- 
posites. This indicates that the hottest dust emission 
is highly extincted or absent in radio galaxi es, as sug- 
gested by orientatio n unification mo dels (e . g..lQgle et al. 
200 1 iClearv et all 120071 lHaas et alJl200a ILeipski et al 
20101) . At A res t < 2 /mi, the radio galaxy SEDs are on 
average about half as luminous as the quasar composite 
for a given A rcs t = 8 /im luminosity density. The radio 
galaxy SEDs at A rcs t < 2 /im are at least as blue as those 
of the quasars and are consistent with a stellar popula- 
tion rather than hot dust emission. In fact, for many 
sources a clear dip is seen between the stellar peak at 
restframe 1.6 /im and the rising dust emission at longer 
wavelengths. It therefore seems feasible to derive robust 
stellar masses from our Spitzer photometry, at least for 
a subset of our sample, though given the large spread 
in SED shapes, great care needs to be taken to separate 
stellar and dust emission in each individual galaxy. The 
next two subsections describe our SED modeling and how 
well we can separate these components. 

3.2. SED modeling 

We used a slightly updated version of the model de- 
scribed by S07 to fit the Spitzer radio galaxy SEDs. For 
the host galaxy component, we used elliptical galaxy 
templates calcul ated from the PEGASEF 5 ! spectr opho- 
tometric model ([Fioc fc Rocca-Volmerangel 11997ft with 
an ass umed formation redshift Zf or m=10 and a iKroupal 
(200.3) initial mass function. (S07 give further details 
about the template and the influences of changing the 
IMF, the formation redshift, and the treatment of ex- 
tinction and thermally pulsing asymptotic giant branch 
stars.) Templates for individual galaxies differ only in 
overall stellar mass and in the amount of passive evolu- 
tion from z = 10 to the individual galaxy redshift. To de- 
scribe the dust emission, we used a simple model consist- 
ing of three pure blackbody components: (i) a T = 60 K 
component, which was used only for the five sources with 
70 /im detections; (ii) a T = 250 K component^; and 
(iii) a component with 500 < T < 1500 K, where the 
temperature was a free parameter. This model with four 
free parameters (five for those sources with 70 /im de- 
tections) was fit to our six (or seven) broad-band data 
points. For sources undetected in the longer wavelength 
data, we fit a maximum dust emission model through 
the 3cr upper limit at A bs = 16 imb. or 24 /im, whichever 
is more sensitive, in order to obtain reliable upper lim- 

13 http://www.iap.fr/pegase 

14 For 6C 0032+412, we had to modify this to a T = 650 K 
component, see §4.3.2. 



its to the hot dust emission. When the two limits were 
within 10 % of each other, the maximum dust emission 
was fitted through both. The contributions from stellar 
populations at A b s = 24 /im are expected to be more 
than an order of magnitude lower than the depth of our 
MIPS photometry and can therefore be ignored. Physical 
dust models with a proper treatment of a clumpy torus, 
viewing angle geom etry, and radiative t r ansfer exist in 
the l it erature (e.g.. iPier fc Krolikl Il992t iGranato et all 
19971 iNenkova et all l2002t lElitzur fe Shlosmanl 120061: 
Fritz et all 120061 : iSiebenmorgen fe Kriigell 120071 ). How- 
ever, a detailed analysis using these models would 
overinterpret our data as the models require mid-IR 
spectroscopy for a d etailed fitting of spectral features. 
iSevmour et all (|2008l ) presented IRS spectroscopy of two 
sources from our sample, and further IRS spectroscopy 
of five sources in our sample will be presented in a fu- 
ture manuscript (J. Rawlings et al. 2011, in prep.). As 
these data cover only 10% of our sample, here we use 
only broad-band photometry and restrict our analysis to 
separating the stellar and hot dust emission in order to 
determine interpolated restframe luminosities. 

Other components besides dust and stars can con- 
tribute to the mid-IR emission. Non-thermal syn- 
chrotron emission has been detected i n Spitzer/IRS spec- 
troscopy of several 3CR sources (e.g JClearv eTaT]l2(M 
ILeipski et all [2009). However, this component is only 
seen in type 1 AGNs dominated by flat radio spectra. 
Our radio galaxies have much steeper radio spectral in- 
dices (see Fig. [TJ) and have low radio c ore dominance (see 
Table [2J|. iMilev fc De Breuckl ((20081 ) showed the SED 
of 4C 23.56, which is one of the most core-dominated 
sources in our sample (CF 2 o = 14.6%), and even for this 
galaxy, the synchrotron contribution near the peak of 
the far-IR SED at A rcs t = 50 /im is at least five orders of 
magnitude lower than the dust emission. We therefore 
do not consider a synchrotron component in our models. 

More important contributions can be expected from 
mid-IR spectral features such as PAH emission at 7.7 /im 
and silicate absorption at 9.7 /im. Both of these features 
have been detected in individual galaxie s in our sample 
([Seymour et al.ll20"08t ILeipski et al.ll2010l ) . However, our 
photometric bands are wide, and the equivalent width 
of PAH emission is low compared to the underlying hot 
dust continuum emission. We therefore do not expect 
these lines to have a significant effect on our photometry. 
The silicate feature may have some discernible effect, but 
as we cover restframe 9.7 /im only for z < 1.5, at most a 
small subset of our sample can be affec ted. Furthermore , 
the composite radio galaxy template of lHaas et all (2008) 
implies that neither PAH emission nor silicate absorption 
is expected to have strong effects on the mid-IR SEDs, 
and thus neither of them is considered in our modeling. 

Finally, the Ha emission line shifts into the IRAC 
3.6 /im band at z > 3.7 (e.g., lAssef et all 120101) . In 
at least one of the seven SHzRG sources at z > 3.7, 
TN J1338— 1942 at z = 4.1, we indeed notice an unusual 
factor of ^2 increase in the 3.6 /im photometry compared 
to the if-band and 4.5 /im points (see Appendix A). We 
therefore use IRAC 3.6 /im photometry at z > 3.7 only 
to verify the consistency of our stellar population but not 
to fit the SED. 

Appendix A shows the SEDs including model fitting of 
the 70 sources in our sample, ordered by increasing red- 
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shift. As already apparent from Fig. [31 the HzRG SEDs 
are quite diverse. Some galaxies such as 7C 1751+6809 
(z = 1.540) show a clear stellar contribution recognized 
by a decline in flux density from restframe 1.6 /xm to 
longer wavelengths and a distinct hot dust component 
steeply rising beyond restframe 3 /im. Other galaxies 
such as 3C 257 (z = 2.474) or PKS 1138-262 (z = 2.156) 
have mid-IR SEDs consistent with a power-law. For 
such sources, we can only derive upper limits to the 
mass of the host galaxy stellar population. The SED of 
TXS J1908+7220 [z = 3.530) is dominated by a trans- 
mitted quasar continuum (see §4.3.2) and also allows only 
an upper limit to be derived on the stellar population. 

3.3. Separation of stellar and dust emission 

Spitzer multi-band photometry is a powerful tool to 
trace the mid-IR emission from AGNs. Various IRAC 
and MIPS color sele ction criteria have been designed to 
identify AG Ns (e. g.. Ilvison et alj|200i lLacv et al.|[200l 
iStern et al.l 12003 iPope et al.l I2008D S07 showed that 
71% and 8 8 % of the SHzRG sou r ces fa ll within the 
IStern et al.l ([20051 ) and lLacv et al.l pOOl IRAC color- 



color AGN selections, respectively. Because 16 sources 
in our sample have only upper limits in the pivotal 8 /im 
b and, the S^^mJS^ ^m versus gs^m/^.sf im diagrams 
of Ilvison et "all (|2004h and lPope etal] (j2008f ) are not very 
useful for our sample. These observed- frame color-color 
selection techniques have two major limitations for our 
purposes: (1) our sample consists, by definition, of type 2 
AGNs, so we are not interested in finding AGNs but 
rather in determining the AGN contribution relative to 
the host galaxy; and (2) our sample covers 1 < z < 5.2, 
while these techniques are not uniformly sen sitive over 
this full redshift range (e.g., lAssef et &U&M \. 

As we have full redshift information for our sample, we 
set up a new restframe color selection criterion to iden- 
tify those sources with strong dust contributions at rest- 
frame 1.6 /im, where the stellar population peaks. To 
interpolate the data between our six photometric data 
points, we used our fitted model (see §3.2). As our bands 
are wide and often contiguous, the use of this model 
should not introduce significant interpolation errors as 
can be seen from the individual SEDs in Appendix A. 
IHaas et all (|2008T ) introduced a restframe <S , 3/zm/>Si.6/im 
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versus S5 pm/Ss urn diagram to separate the quasar and 
radio galaxy SEDs in their 3CR sample at 1 < z < 2.5. 
We modified this criterion for use at z > 2.5, where MIPS 
24 pm photometry only covers out to A rcs t ^$ 5 /im. On 
the blue end, the data cover the peak of the stellar emis- 
sion at A rcs t < 1.6 /im for our entire sample. Consider- 
ing an even bluer A rcs t would be possible observationally 
but would increase the contributions from younger stellar 
populations (see S07) as well as non-stellar contributions 
suc h as scattered quasar emi ssion and nebular emission 
(see lMilev fc De Breucldl2008j ). We also wanted our crite- 
rion to pick up a possible minimum in the SEDs longward 
of the Rayleigh- Jeans fall-off of the stellar emission at 
Arcst > 2 pm and the onset of the steeply rising, hot dust 
continuum. While the stellar peak has a unique shape, 
the dust continuum slope depends on the temperature 
and obscuration of the hot dust emission. Moreover, the 
relative offset between the stellar and hot dust compo- 
nents can vary substantially as discussed in §3.1. We 
therefore identified A ros t = 3 pm as the pivot wavelength 
for our new criterion. This wavelength is also located 
centrally between 1.6 and 5 pm, which minimizes the ef- 
fect of measurement errors on our photometry. Table |4] 
lists the monochromatic luminosities at 1.6, 3, and 5 pm 
as derived from our SED fitting. Table 0] also lists the 
fraction of stellar light at A re st =1.6 /im and the stellar 
masses derived following the procedures of S07. 

Figure 0] shows the restframe 5 , 3^ tm /S'i.6/xm ver- 
sus S$ S3 am color-color diagram for our sam- 
ple. All AGN-dominated SEDs cluster around 
(<S , 5/im/S , 3 ^ m , -S , 3 /ffli/S'i. 6 /mi) ~ (1.5,3). This corre- 
sponds to an SED which is still steepening at A rcs t < 
3 /im, indicative of the warm temperature limit of the 
dust emis sion with a possible contribution fro m extinc- 
tion (e.g. JHaas et al.ll2008t ILeipski et al.M2010f ). Sources 
with very red SEDs at A res t < 3 /im (upper left corner in 
Fig. [4J are therefore likely to be dominated by AGN dust 
emission. Pa rk et all (|2010l ) showed that Spitzer power 
law galaxies have 3-24 /im spectral indices — 1.8 < a < 
(S v oc v a ). We therefore defined the straight power law 
locus as a boundary in Fig.0]to separate AGN dust emis- 
sion and host galaxy (stellar) dominated HzRGs. Sources 
with a < —1.8 are indicative of strong extinction already 
in the A re st = 3-5 /im range and are not expected to have 
significant remaining dust emission at 1.6 /im. Sources 
with S3 urn/ S\.q < 1 are consistent with stellar pop- 
ulations as confirmed by the positions of the PEGASE 
and M82 SED in Fig. 01 Summarized, we defined sources 
dominated by an AGN at 1.6 /im by the criterion: 

S3 pm/ S\.q pm > 1 

n 

S3 /im/ Si. & fim > 1-39 x S5 p,m/ S3 /im — 0.39 

n 

S5 /im/ S3 fjm < 2.5 

Even with our criterion, classification of individual ob- 
jects is uncertain, especially near the boundaries: (1) 
measurement errors may shift objects across the bound- 
aries, (2) variations in the interpolation between the 
bands using the SED model described in §3.2 cause 



some uncertainties, and (3) contributions from other non- 
stellar components such as a transmitted or scattered 
quasar continuum might require a more complex model 
for some sources (e.g., TXS J1908+7220, see §4.3.2). For 
individual sources, a more accurate separation would be 
possible when using additional supporting data such as 
IRS spectra or multi-band near-IR and visible imaging. 
This is beyond the scope of the current paper and would 
also introduce redshift dependence in the results. We 
are mainly interested in looking for statistical correla- 
tions, and our rest-frame color criterion allows a uniform 
treatment for all sources within our sample. 

For 15 sources in our sample (mostly at z > 2), we 
could obtain only upper limits on the hot dust emis- 
sion; these were determined from the MIPS 24 pm limits 
only, as our MIPS images are significantly more sensi- 
tive than the IRS 16 /im imaging. However, the IRS and 
MIPS photometry is sufficiently sensitive to assure that 
the IRAC detections are dominated by stellar emission 
(e.g., TN J1338-1942 at z = 4.1). In total, we derived 
reliable stellar masses for 42 HzRGs (see also §4.1). For 
the remaining 28 sources, we used the same procedure to 
correct the A rcs t = 1.6 /im flux from the extrapolated hot 
dust model and fit the stellar template to the remaining 
photometry. This procedure does become less reliable: 
the median 1.6 /im stellar fraction of these 28 sources is 
0.71 compared to 0.96 for the sources outside the AGN- 
dominated area (see tabled]). While the stellar masses 
may be close to the real ones (especially for sources near 
the boundary in Fig. [4]), we conservatively treat them as 
upper limits. 

4. RESULTS 

Having identified that both stellar mass and AGN- 
heated dust emission contribute significantly to the near- 
and mid-IR SEDs of HzRGs, we now examine the cor- 
relation between the physical parameters in our sample. 
In order to properly consider the upper limits in the sta- 
tistical tests of the significance of correlations, we make 
use of survival analysis and i n particular the ge neralized 
Spearman's p coefficient (e.g.. Hsobe et al.|[l986D . Results 
are expressed in terms of the probability P of seeing the 
observed or a higher value of p if the null hypothesis that 
the data are uncorrelated holds. Thus a small value of 
P indicates a high likelihood that a correlation exists. 

4.1. Stellar masses 

The original motivation of our Spitzer survey was to 
derive accurate stellar masses by consistently observing 
the same restframe wavelength where the old stellar pop- 
ulation peaks, which S07 argued is at restframe .ff-band. 
Figure [5] updates the HzRG restframe _ff-band stellar lu- 
minosity density-redshift relation. The only difference 
from the S07 version is that with full IRS and MIPS 
photometry, we now have a better handle on the hot 
dust emission, particularly for sources where previously 
we obtained only upper limits on the stellar masses (see 
§3.3). The stellar rest-frame iJ-band luminosity shows 
no correlation with redshift (P = 0.42). 

Figure [5] shows the redshift evolution of the derived 
stellar masses. We confirm the result of S07 that HzRGs 
have high stellar masses at every redshift. However, we 
now find a trend for the z > 3 radio galaxies to be slightly 
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Figure 5. Rcstframe H-band stellar luminosity versus redshift for the Spitzer HzRG sample, derived from the best-fit models to the 
multi-band photometry. Upper limits indicate radio galaxi es where only a maximum fit to the stellar SED was possible. The dashed lines 
show the luminosities of PEGASE.2 elliptical galaxy models ( Rocca-Volmcrangc ct al. 2004) with -Sform — 10 and masses 10 11 and 1O 12 M . 



less massive (P = 0.0044). The drop in stellar mass be- 
comes apparent only at z > 3. Considering only the 
subset of sources at z < 3, the correlation becomes in- 
significant (P = 0.19). We caution that the putative 
mass decrease is based on only 11 reliable stellar mass 
determinations at z > 3, and the decrease in stellar mass 
is very small: the median stellar mass of z < 3 HzRGs 
is 2.3 x 10 11 M versus 1.6 x 10 11 M for z > 3. If con- 
firmed, the difference would indicate that radio galaxies 
are still forming at z > 3 but have already assembled the 
bulk of their stellar mass by z ~ 3. This pivotal red- 
shift is consistent with the increase in the submillimeter 
detection rate of radio galaxies at z > 3, where submm 
flux densities above a few mJy imply e xtreme star for- 
mation rates exceeding 10 00 Mq yr" 1 (| Archibald et al.l 
[20011 IReuland et all 120041. A more elaborate spectral 
decomposition would be required to better quantify the 
reality of this small drop in stellar mass. In particu- 
lar, rest-frame visible and UV photometry and polarime- 
try are needed to properly extrapolate the scattered and 
transmitted AGN contributions. In addition, dust ob- 
scuration may affect even restframe n ear-IR emission in 
the host galaxy (e.g. JDev etalll2008f) . While the bright 
submm emission in some z > 3 radio galaxies suggests 
they may indeed contain large amounts of dust (though 



the submm brightness is most likely dominated by higher 
star formation rates), we do not expect this dust to 
cause large obscurations because HzRG extended emis- 
sion lines, in p articular Lyq, are less ob scured at z > 3 
than at z < 3 dVillar-Martm et alj|2007ft . 

4.1.1. Dependence of stellar mass on radio luminosity and 
black hole activity 

Figure [7] revisits the marginal correlation between stel- 
lar mass and radio power discussed by S07. Our new 
plot is more reliable in two ways: (i) it uses the more 
nearly isotropic 500 MHz rather than the 3 GHz radio lu- 
minosity, and (ii) it has more reliable determinations of 
stellar mass upper limits (see §3.3). Our full data set 
confirms the result of S07 that radio power and stel- 
lar mass are not significantly correlated (P — 0.26). 
The absence of a stronger M s t e iiar _ i500 MHz correlation is 
slightly surprising given previous claims of a correlation 
between radio power and observed if -band magnitude 
(lEales et al.l Il997t IBest et all [l99l IJarvis et al.l l200ll: 
IDe Breuck et al.ll2002t Twillott e t alll2003D . Apart from 
the difficulty in disentangling redshift and radio luminos- 
ity effects in samples affected by a strong Malmquist bias, 
the correlations with if -band magnitude apparently stem 
from other components that correlate with radio power 
such as transmitted and/or scattered quasar light, neb- 
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Figure 6. Stellar mass versus redshift for the Spitzer HzRG sam- 
ple. Masses are based on the SEP decomposition describ ed in >B1 
and the PEGASE.2 models IIRocca-Volm crangc ct al. 200?). 



ular emission (continuum and/or lines), or (in the low- 
redshift objects) the tail of the hot dust emission. 
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Figure 7. Stellar mass (ij4} versus restframe 500 MHz radio lu- 
minosity for the Spitzer HzRG sample. 

The limited radio luminosity range of our sample and 
the small spread in stellar masses make it difficult to 
detect any but a strong correlation. However, we can 
also turn this argument around and suggest that the 
low stellar mass spread is proof that radio power, stel- 
lar mass, or both are not strongly correlated with the 
accretion properties of the central supermassive black 
hole (SMBH). The correlation between the mass of the 
nuclear bl ack hole and the stellar bulge is now well es- 
tablished (|Magorrian et al.lll99a fearing fc Rixll2004D . If 
this relation also holds for HzRGs, the small range in 
stellar masses would also imply a small range in black 
hole masses. Based on spectroscopy of six broad Ha 



lines in HzRGiO, iNesvadba etal] (|2010l) do indeed find 
all black hole masses in the range A/bh = 3 x 10 9 to 
2 x 10 10 M©, which puts them only slightly offset from the 
local MsH-Aibuige relation. While further Mbh determi- 
nations would be needed to check this trend, it seems 
likely that HzRGs contain SMBHs with masses similar 
to those expected. 

On the other hand, low frequency radio luminos- 
ity is only one of the e s timators of the jet k inetic 
power (jWillott et al.l [19991 ICattaneo fc BestH2009D and 
hence black hole activity. Comparisons with other 
measurements, such as the amount of kinetic power 



needed to create X-ray cavities (jBirzan et al.H2008l I2010t 
ICavagnolo et al.|[2010j) . show that there can be an un- 
certainty of more than an order of magnitude in this 
determination. Intrinsic spread in the radio luminosity- 
jet power relation, environmental effects (e.g., density of 
the medium in which the radio source is expanding), and 
variations in the Eddington ratios of SMBHs of similar 
mass may therefore imply that the 500 MHz luminosity 
does not provide a very reliable measure of the black hole 
properties, at least not in the limited 2 dex range of the 
most luminous HzRGs studied by our sample. 

4.2. Mid-IR luminosities 

Before the launch of Spitzer, AGN-heated dust 
emission remained a virtually unexplored component 
in z > 1 radio galaxies, apart from mostly up- 
per li mits obtained by the Infr a red Space Observatory 
(e.g.. iMeisenheimer et al.l 120011 ; lAndreani et al.l 120021; 
lHaas et al.l I2004D . S07 presented Spitzer results regard- 
ing the mid-IR luminosities of HzRGs, and Section 3 ad- 
vances the argument that the A rcs t > 3 /im emission of 
HzRGs is dominated by hot dust. The complete IRS 
16 /im and MIPS 24 /im photometry presented herein 
permit further analysis. Hot dust emission was detected 
in 57 out of 70 HzRGs; the non-detections are mainly 
at z > 3 and likely are a result of using a fixed integra- 
tion time per object. This section examines the charac- 
teristics of the hot dust emission and what it can reveal 
about AGN activity. We selected a restframe wavelength 
of 5 /im in order to use a consistent luminosity through- 
out the redshift range in our sample. Contributions from 
stellar emission are negligible at this wavelength (see Ap- 
pendix A). 

4.2.1. Orientation effects 

Low-frequency radio emission is considered to be one 
of the most nearly isotropic components seen in distant 
AGNs. In contrast, mid-IR emission originates from the 
innermost regions of the torus and is expected to be much 
less isotropic due to obscuration effects. The mid-IR 
luminosity in type 2 AGNs does indeed appear to be 
fainter (= more obscured) than in type 1 AGNs from 
the same parent radio-selected samples, provid ing strong 
support for the orientation unification m odel (Sh i et al.l 
2005; lHaas et al.ll200l ILeipski et al1l201Cl) . 

While our sample is also selected at low frequency, it 
contains only type 2 AGNs as determined from their 
restframe ultraviolet (observed visible) spectra. How- 
ever, our mid-IR SEDs show a large variety, ranging 

15 Five of which are also in our SHzRG sample. 
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from clearly host-galaxy-dominated sources to SEDs that 
are almost indistinguishable from those of type 1 AGNs 
(Fig. [3]). This suggests a gradual transition from highly 
obscured type 2 AGNs all the way to unobscured type 1 
AGNs. To verify this hypothesis, we have checked 
whether radio morphological parameters can be used 
as indepe ndent indicat ors for orientation. Both radio 
size (e.g., [B arthcl 1989) and radio core dominance (e.g., 
lOrr &: Brownelll982D have been suggested as measures of 
orientation. The largest radio size, 6, and the hot-dust- 
dominated L^^ui show a likely correlation (P = 0.014). 
This suggests that 9 is probably an indicator of orien- 
tation but not a very reliable one. On the other hand, 
Figure [8] shows that the core fraction at restframe fre- 
quency 20 GHz, CF20, is strongly correlated with Ls^trn 
(P < 0.0001). As we could determine core fractions of 
only 34 of our 70 sources, we caution that this analy- 
sis is not complete. Most of the 25 double-component 
sources have radio cores too faint to be detected in 
the current radio maps and are expected to have core 
fractions CF20 << !%• These 22 sources (marked by 
a 'D' on the righthand side of Fig [8]) have a median 
(log(L5 iim/Lo)) = 11.7 and will fall close to the existing 
correlation. The remaining 8 single-component sources 
have a median (log(L 5 ^ m /L )) = 11.8, but without 
higher spatial resolution observations we cannot estimate 
their Ci^o- However, they represent only 11% of the 
sample and are unlikely to be able to remove the correla- 
tion completely. We now examine the physical origin of 
this correlation between the synchrotron-dominated ra- 
dio emission and the hot dust dominated A rcs t = 5 jum 
emission. (As argued in §3.2, a direct contribution from 
synchrotron emission at A rcst = 5 fim can be neglected.) 



13 



12 



10 



++ 

t + + + 

I + + 

+ + 



+ + J 
I + 



rect view of the inner parts of an obscuring torus perpen- 
dicular to th e radio jet. This inte rpretation was recently 
confirmed bv lLeipski et al.l ((2010), who used Spitzer/IRS 
spectroscopy to show that the silicate absorption depth 
79.7 /xm increases with decreasing radio core fraction (but 
see Landt et al. 2010 for a diff erent view). Ba sed on a 
sample of 31 3C radio galaxies. lOgle et all (|2006[ ) claimed 
that radio morphology is not a reliable predictor of nu- 
clear mid-IR luminosity. However, considering only their 
13 "mid-IR luminous" radio galaxies, there does seem 
to be a correlation, consistent with our results which 
are based on a sample three times larger. If the ra- 
dio jets in the more core-dominated and more luminous 
hot dust sources are indeed oriented closer to the line 
of sight, we would expect other anisotropic parameters 
to increase also. Examples include the equivalent widths 
of the [O II] A3727A and [O III] A5007A lines, both 
of which are known to correlate with radio core domi- 
nance in type 1 AGNs (jBaker fe Hunsteadlll995D . Deep 
restframe visible spectroscopy of our sample is currently 
underway and should be able to test this hypothesis (Nes- 
vadba et al. 2011, in prep.). 

We conclude that the radio core fraction is a powerful 
indicator of the expected mid-IR luminosity because it 
provides a measure of the obscuration. As a result, the 
spread in L5 um in our sample is in large part driven 
by obscuration/orientation effects with the most core- 
dominated sources being the least obscured ones. 
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Figure 9. Monochromatic luminosity at 5 fim versus redshift. 
Horizontal bars indicate median luminosities in each of three red- 
shift bins. Upper limits, indicated by arrows, cluster at z > 3 
owing to the use of a fixed integration time for all sources. 



Figure 8. 5 (im luminosity versus radio core fraction at restframe 
frequency 20 GHz. The 5 fim luminosities of the double- and single- 
component radio sources that lack measurements of their core frac- 
tions are indicated on the right hand side of the plot by D and 
S respectively. The evident correlation indicates an orientation- 
dependence in our sample. 

Higher core fraction sources are more strongly beamed, 
which can be understood if the radio jet is oriented closer 
along the line of sight of the observer (e.g. JOrr fc Browne! 
1 19821 ). These sources are therefore more similar to type 1 
AGNs. The higher hot dust luminosity in them (Fig. [8} 
also indicates less obscuration, consistent with a more di- 



4.2.2. Evolutionary effects 

Figure M suggests a trend of rising restframe 5 /im lumi- 
nosity with redshift, which is confirmed by our general- 
ized Spearman rank analysis (P = 0.0021). Before inter- 
preting this correlation in terms of evolution of the dust 
content or emissivity, it is important to consider three 
selection effects. First, the 5 /im luminosity is even more 
strongly correlated with radio power (see below), and 
there is a small remaining Malmquist bias in our radio- 
selected sample (see §2.1). This combined effect can am- 
plify the redshift evolution of L5 though it is unlikely 
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to imprint a correlation as strong as observed. Second, as 
shown in §4.3.1, orientation effects can strongly influence 
L^fim- However, the core fractions are distributed quite 
evenly in redshift, suggesting that there is no redshift- 
dependent obscuration effect. Third, 9 out of 11 upper 
limits on L5 are at z > 3 compared to 13 detections 
at z > 3. As shown by the survival analysis, the data 
still indicate a real correlation, though deeper mid-IR 
photometry at z > 3 will be needed to quantify the ac- 
tual strength of the correlation. If confirmed, this would 
indicate that higher- redshift AGNs are more efficient in 
heating dust. 

4.2.3. Correlation with radio power 

Several authors have compared the hot dust emission 
from distant radio-loud AGNs at 8 < A rBS t. < 30 /nm 
with the 178 MH z radio luminosity (e g..lShi et al.ll2005b 
lOgle et al.1 120061 iClearv et all 120071: lHaas et all 120081 ). 
They found that type 1 AGNs show a clear positive 
correlation with a scatter of roughly an order of mag- 
nitude. Type 2 AGNs show more scatter and generally 
fainter hot dust emission because of obscuration, though 
the upper end of the type 2 range reaches the locus 
of the type 1 AGNs. Figure \W\ plots 5 /Ltm luminosity 
against 500 MHz radio luminosity density for our sam- 
ple and shows a strong correlation (P < 0.0001). The 
correlation is equally significant when using the 3 GHz 
luminosity instead. Obscuration will scatter sources 
downwards in Fig. 1101 and it is remarkable that the 
observed correlation is so strong despite the orienta- 
tion effects discussed in §4.3.1. We can get an idea 
of the intrinsic scatter by considering only the sources 
with logL 5 ^ui/Lq > 12.5, which are likely to experi- 
ence only minimal extinction. These 8 sources cover a 
range of radio powers 28.3 < logLsooMHz < 29.3. If 
all sources within this radio power range have the same 
intrinsic logLs^m, they exhibit up to 2 orders of magni- 
tude of extinction. The fact that none of the 23 sources 
with logL 500 MHz < 28.3 has log L 5 ^ m /L Q > 12.5 indi- 
cates that the ratio L5 ^m/ £500 MHz has an upper limit 
around a few hundred (see dotted lines in Fig. H0| . How- 
ever, it is impossible to determine the intrinsic spread 
in £5^ul/L5oomhz because of the range of extinction in 
Ls^m- Orientation effects will increase the spread in the 
i5/im~^500MHz plane, and the intrinsic correlation may 
be even stronger than observed. 

The L5 /j I m-L5Qo MHz correlation implies that both are 
related to the energy output of the SMBH despite the 
vastly different spatial scales of the emission. The hot 
dust emission originates from the inner regions of the 
host galaxy near the S MBH on scales w hich may be as 
small as a few pc (e.g.. Uaffe et al.l[2004[ ). while the radio 
luminosity is dominated by the extended radio lobes of- 
ten extending beyond the host galaxy and reaching scales 
up to hundreds of kpc (Tabled]). It is therefore remark- 
able that the scatter in the relation is limited to only 
two orders of magnitude. Naively, one might have ex- 
pected a much larger spread, as it may be influenced by 
several effects such as (i) orientation effects as discussed 
above, (ii) different Eddington accretion rates onto the 
SMBH, (iii) time delays related to the very different size 
scales of the emitting regions, (iv) different dust distribu- 
tions or compositions (e.g., in 6C 0032+412, see §4.1.2), 
and (v) environmental effects related to the medium in 
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Figure 10. Restframe 5 fim luminosity against monochromatic 
500 MHz luminosity. Dotted lines indicate mid-IR to radio ratios 
of 1, 10, 100 and 1000 (bottom to top). 

which the extended radio sources are expanding. The 
fact that the radio luminosity is more strongly corre- 
lated with hot dust emission than with the stellar mass 
of the host galaxy (§4.2.1) suggests that variations in the 
Eddington ratio of the AGN may be the dominant fac- 
tor in driving both t he ra dio p ower and hot dus t emis- 
sion. iHickox et all (f2009T) and iGriffith fc Sterol (|2010f ) 
also suggest that radio-loud AGNs show a larger range 
in Eddington luminosity than X-ray and mid-IR selected 
AGNs. 

4.2.4. Absence of a correlation with stellar mass 

Does hot dust emission depend on the mass of the host 
galaxy? Figure [TT] plots the hot dust dominated 5 /xm lu- 
minosity against stellar mass, showing that the two are 
uncorrelated (P = 0.73). The upper limits in the top 
right corner illustrate the success of our color-color selec- 
tion criterion to identify hot dust dominated SEDs (§3.3, 
Fig. [4|. The absence of a correlation also suggests that 
AGN luminosity is independent of stellar mass, though 
we warn that orientation effects (§4.2.1) make is^tm not 
a reliable measure of AGN luminosity. 

4.3. Notes on individual sources 

Radio galaxies are often complex systems consisting 
of several components. Therefore, our Spitzer photome- 
try is occasionally affected by companion and/or fore- 
ground objects. We have visually inspected the six- 
band Spitzer images for our entire sample, using multi- 
wavelength data in the literature to determine the correct 
astrometric identifications as listed in Table [TJ In most 
cases, we found an unambiguous radio source counter- 
part. The next two subsections describe all special cases, 
first sources with multiple AGNs (§4.3.1) followed by de- 
tailed discussions of a few additional interesting sources 
(§4.3.2). 

4.3.1. Binary AGNs 
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Figure 11. Stellar mass versus 5 fim luminosity. The lack of 
correlation implies that our criteria to isolate hot-dust dominated 
SEDs are successful and that stellar mass is uncorrelated with 5 (im 
luminosity. 

As shown in §3) HzRGs are massive galaxies hosting 
a SMBH. If such massive galaxies are formed in a series 
of merger events as predicted by the standard A Cold 
Dark Matter paradigm, we expect that merging massive 
galaxies will coincide w ith merging central SMBHs (e.g., 
IComerford et al. 2009). Observationally, this may man- 
ifest itself as a second AGN near the HzRG, or — at 
later stages in the merger event — by perturbed mor- 
phologies as observed in low redshift radio galaxies (e.g., 
iRamos Almeida et al.l[20 10). The low spatial resolution 
of Spitzer, particularly at the longest wavelengths, limits 
such studies to companion AGNs at many-kpc distances 
rather than the sub-kpc separations in later stages of 
merging. The mid-IR images show evidence for a second 
AGN in four HzRGs in our sample (Fig. [T2] discussed in 
RA order): 

4C 60.07 (z = 3.788): llvison et all d200l showed 
that this system consists of two AGNs separated by 4" 
(30 kpc) . The radio galaxy (A) , as identified by the radio 
core, remains undetected in the hot dust emission beyond 
4.5 /im, while the companion object B, which currently 
lacks a spectroscopic redshift, has very red hot dust emis- 
sion characteristic of an obscured AGN. An initial claim 
of a C O(l-O) redshift for this second AGN (|Greve et al.l 
120041 ) could not be conf irmed in the re-analysis of those 
data (jlvison et al J 12008). but the (sub) mm dust contin- 
uum morphology looks consistent with both targets be- 
ing at the same redshift. 

3C 356 (z = 1.079J: There are two potential identifi- 
cations for the location of the radio-loud AGN. Based on 
the detection of scattered broad line regions with a polar- 
ization ang le perpendicular to the "dum bell" structure of 
the galaxy, iSimpson fc R awlingsj (|2002f) argued that the 
northwestern object a is the most likely location, but ob- 
ject b, offset by 4'.'5 (36 kpc), may also contain an AGN. 
Both target s were spectroscopica lly confirmed to be at 
z = 1.079. iCimatti et al.l (|1997| ) showed that object a 
has a high restframe UV polarization, clearly indicating 
the presence of a hidden AGN. We detected both objects 
in all four IRAC bands, though only object a has strong 
mid-IR dust emission seen with IRS and MIPS. This pro- 
vides further support for the conclusion that object a is 



the location of the radio galaxy host, while object 6 is a 
companion object. Tables [3] and [5] list the photometry 
for both objects. 

7C 1756+ 6520 (z = 1.416J: IGalametz etafl (pOOflh 
found three objects within 6" (48 kpc) of the radio 
galaxy, two of which w ere spectroscopically confirmed by 
IGalametz et al.l (|2010h to be at the same redshift as the 
radio galaxy. One of them is a mid-IR selected AGN. 
Both this AGN and a passively evolving BzK-selected 
(pBzK) galaxy were individually detected with IRAC 
but only marginally resolved with IRS and MIPS. Ta- 
bles [3] and [5] provide photometry of both objects; as the 
HzRG and AGN are not fully spatially resolved in the 
IRS 16 fim and MIPS 24 /im images, we estimated the 
relative fluxes between the HzRG and AGN by summing 
the pixel values in a fixed circular 6 pixel radius aper- 
ture using the positions from the IRAC 3.6 /im image 
(Fig. [12]) as centroids. The derived ratios of 50%/50% at 
16 /im and and 20%/80% at 24 /im should be considered 
approximative. This paper use s the r evised redshift of 
the HzRG from IGalametz et al.l (|2010[) . 

MRC 2048-272 (z = 2.060J: No radio or X-ray core 
was detected in thi s radio galaxy (iCarilli et al] 119971 : 
lOverzier etHI 120051 ). iPentericci et al.l (|2001[ Tdetected 
two sources, separated by 2'.'5 (21 kpc), in NICMOS H- 
band imaging and considered the easternmost one as 
the most likely AGN host. This is consistent with the 
Lycn halo, which exten ds for ~5" along the radio axis 
(jVenemans et al.l 120071 ). We detected both components 
in the IRAC 3.6 and 4.5 /xm images. However, at longer 
wavelengths only the western source is detected, suggest- 
ing it also contains an AGN. This configuration is remi- 
niscent of 4C 60.07. 

4.3.2. Other Sources 

6C 0032+412 (z = 3.670): This HzRG is the only 
one where our standard 3-component dust model (using 
T = 60, 250, and 500-1500 K, see §3.2) could not fit the 
SED. First, we had to increase the temperature of the 
T = 250 component to T = 650 K. In addition, this is the 
only source where we required T = 1500 K for the hottest 
dust component compared to a median (T) = 800 K . A 
near-IR bump caused by a T ~ 1500 K dust component 
has been seen in several quasars (e.g.. iGallagher et al.l 
120071: iMor et aJJl2009t iLeipski et al.ll2010t) and is gener 
ally associated wit h graphite dust near the sublimation 
temperature (e.g.. IBarvainislll987| ). The very high tem- 
perat ure necessarily locat es the dust very close to the 
AGN. lLeipski et aLl (|2010f ) argued that this very hot dust 
emission is either more highly extincted or completely ab- 
sent in radio galaxies. While this object is the only ex- 
ample in our sample, it does represent the first counter- 
example to the radio galaxy trend. 6C 0032+412 is a 
bona-fide type 2 AGN as derived from its nar row-line UV 
and visible spectroscopy (jJarvis et al.ll200lT ). The mod- 
est 2% radio core fraction (Tabled also does not indicate 
a line of sight close to that of quasars (see §4.3.1). A pos- 
sible explanation would be that 6C 0032+412 presents a 
rare unobscured line of sight through the entire dusty 
torus all the way to the central regions near the AGN. 

MRC 0037-258 (z = 1.10 j: The new radio imaging 
clearly detected the radio core at 4.86 and 8.46 GHz. 

MRC 0114-211 (z = 1.41 j: We identify the faint cen- 
tral component, which is spatially resolved only in the 
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Figure 12. Spitzer /IR AC 3.6 /jm images of the four binary AGNs in our sample. Contours show the MIPS 24 /j,m images. 
Photometry for the HzRG and companions are listed in Tables [3] and O respectively. 



8.46 GHz map, as the core. 

6C* 0132+330 (z = 1.71): We identify the central 
component seen at both freq uencies as the radio cor e. 

6C 0140+326 (z = 4.41 j: IRawlings et all (fl996h ar- 
gued that this radio galaxy is gravitationally lensed by a 
foreground (z = 0.927) galaxy 1'.'6 to the southeast. This 
foreground L» galaxy has an e llipsoidal shape in both 
the Keck/NIRC if -band image (Ivan Breueel et allll998l ) 
and the i/ST/NICMOS F 1 60 W image (lLacvHl999lh Our 
Spitzer photometry at 3.6 and 4.5 /jm is dominated by 
strong a strong foreground star, while the 16 and 24 /im 
photometry are dominated by mid-IR emission from the 
lensing galaxy. The radio galaxy itself appears unde- 
tected in all but the 5.8 or 8 /im emission, where it is 
only marginally detected and contaminated by the fore- 
ground galaxy. We therefore provide only relatively shal- 
low upper limits based on the detection of the foreground 
galaxy. 

MRC 0251-273 (z = 3.16): Our radio map (Ap- 
pendix B) shows five bright components at 8.46 GHz 
as well as diffuse emission to the north-northwest 
(mainly at 4.85 GHz ). The lower- resolution radio map of 
iKapahi et al.l (|1998t ) confirms that this diffuse emission 
is real. Within the radio structure, there is only a sin- 
gle IRAC identification, coinciding within 0'.'3 with the 
northernmost of the bright components in the 8.46 GHz 



map. We therefore identify this component as the ra- 
dio core. T his radio m orphology resembles that of 
B2 0902+34 (|Carillil H995I ). The radio jet is probably 
aligned relatively closely to the line of sight, with the 
complex south-southeastern radio structure being the ap- 
proaching radio jet. 

MRC 0324-228 ( z = 1.894): The NIC MOS image 
of this radio galaxy (jPentericci et al.ll2001| ) reveals two 
potential host galaxies separated by 1'.'5 (12kpc). Our 
Spitzer imaging does not resolve them as individual 
sources, though it does show an extended structure along 
the same north-south orientation. Our new radio data 
tentatively detected a faint radio core at 4.86 GHz, co- 
inciding with the southernmost of these two identifica- 
tions, which we have adopted as the host galaxy. There 
are also two nearby objects that may contaminate the 
MIPS photometry. 

MRC 0350- 279 (z = 1.90J: No radio core is detected 
in our new radio maps. 

5C 7.269 (z = 2. 218/. No radio core is detected in our 
new radio maps. 

6C 0820+3642 (z = 1.860J: The new radio imaging 
detects a faint radio core, which is most obvious in the 
8.46 GHz map. 

6C 0901+3551 (z = 1.91 No radio core is detected 
in our new radio maps. 
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USS 1243+036 (z = 3.570): Our Spitzer photometry 
is not sufficiently sensitive to exclude a possible contri- 
bution from hot dust at rest-frame 1.6 /im. The derived 
stellar masses are therefore formally considered as upper 
limits, although the 3.5 and 4.5 /im photometry suggests 
they are most likely uncontaminated by hot dust. 

USS 1707+105 (z = 2.349;.- The new radio imaging 
detects a faint radio core, which is most obvious in the 
4.86 GHz map. 

LBDS 53W002 (z = 2.293;.- The IRS and MIPS pho- 
tometry are po ssibly contamina ted by a galaxy ^4" to 
the northwest (|Keel et al J 12002). This galaxy is most 
likely at a different rcdshift than the AGN host as it does 
n ot appear in the [ O III] and Ha narrow-band imaging 
of lKeel et al.l p002) . 

7C 1751+6809 (z = 1.54): No radio core is detected 
in our new radio maps. 

7C 1805+6332 (z = IM): The new radio imaging de- 
tects a faint radio core, which is only seen in the 8.46 GHz 
map. 

TXS .1190 8+7220 (z = 3.530 J: Using integral field 
spectroscopy, iSmith et~a l. (2010) argued that this broad 
line radio galaxy is a system of two vigorously star- 
forming galaxies superimposed along the line of sight and 
separated by ^130 kms" 1 in v e locity . The detailed SED 
decomposition of ISmith et al.l ((20100 also showed that 
the ncar-IR emission of this HzRG is dominated by trans- 
mitted quasar continuum. While the hot dust continuum 
does not contribute substantially at A rcs t = 1-6 /im, this 
transmitted quasar continuum prevents us from measur- 
ing the old stellar population, and we provide only an 
upper limit to the stellar mass. 

TN J 2007- 13 16 (z = 3.84;.- We identify the second 
component from the north as the radio core. This com- 
ponent is spatially resolved only in the 8.46 GHz map. 

MG 2144+1928 (z = 3.592;.- Our Spitzer photometry 
is not sufficiently sensitive to exclude a possible contri- 
bution from hot dust at rest-frame 1.6 /mi. The derived 
stellar masses are therefore formally considered as upper 
limits, although the 3.5 and 4.5 /im photometry suggests 
they are most likely uncontaminated by hot dust. 

5. SUMMARY AND CONCLUDING REMARKS 

The main results from our full six-band Spitzer imag- 
ing survey of 70 radio galaxies at 1 < z < 5.2 can be 
summarized as follows: 

1. Using a newly designed criterion to isolate sources 
dominated by hot dust emission at A rcs t=1.6 /im, 
we unambiguously detect a stellar population in 46 
radio galaxies and hot dust emission in 59. 

2. Four of the 70 HzRGs show a second AGN within 
6", as revealed by their mid-IR colors. This may 
be an indication of an ongoing or imminent major 
merger event. 

3. The stellar masses have a remarkably low scatter 
around ^2 x 10 11 M Q . There is tentative evidence 



for a small drop in stellar mass at z > 3. If con- 
firmed by more detailed spectral decompositions, 
this result suggests that radio galaxies have built 
the bulk of their stellar populations by z ~ 3. This 
is consistent with the observed increase in their 
submillimeter luminosities, suggesting higher star 
formation rates, at z > 3. 

4. The radio core fraction is strongly correlated 
with restframe 5 /tm luminosity. The most core- 
dominated radio galaxies are also the most lumi- 
nous (i.e., least obscured) hot dust emitters. This 
suggests a gradual increase of the hot dust obscura- 
tion from type 2 to type 1 AGNs, consistent with 
the orientation paradigm of AGN unification. It 
also complicates the use of A rcs t = 5 /im luminosity 
as an indicator of AGN power. 

5. The 500 MHz radio luminosity is only marginally 
correlated with stellar mass but strongly correlated 
with the hot dust luminosity at A rcs t = 5 /im. Vari- 
ations in the Eddington ratio of the SMBH could 
explain these trends. 

Our Spitzer sample provides an unique database to 
study the massive host galaxies of powerful radio galax- 
ies, their AGN-heated dust emission, as well as their envi- 
ronments. In future papers we plan to study the stellar 
populations in greater detail by combining the Spitzer 
data with visible and near-infrared data to better con- 
strain the HzRG star formation histories and thus derive 
more accurate stellar masses and ages. We also plan 
to observe the cold dust component with the Atacama 
Large Millimeter and submillimeter Array (ALMA) and 
the peak of the dust emission with Herschel Space Ob- 
servatory. This full SED will allow to separate the hot 
AGN from the cool, starburst-heated dust emission, pro- 
viding more accurate measures of the instantaneous star 
formation rates. This combined dataset will provide a 
unique view of the star formation history in these very 
massive galaxies at every epoch. 
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APPENDIX 

SED FITS TO RADIO GALAXIES AND NOTES ON INDIVIDUAL SOURCES 

Figure [13] presents the Spitzer SEDs for our complete HzRG sample, including model fits to the mid-IR photometry, 
ordered by redshift. 
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Figure 13. Radio galaxy SEDs and model fits. Lower abscissas arc marked with observed wavelengths, and the scales arc the same for all objects. 
Upper abscissas arc marked with rcstframc wavelengths, which depend on wavelength for each object. Filled circles with error bars denote data 
points used in the fit; open circles denote data points that could be contaminated by emission lines and arc therefore not used in the fits. Model 
components are a stellar population (green dashed lines) and two or three pure blackbodics (red dotted lines) representing dust emission. The sum 
of the model components is shown by solid black line. The name and rcdshift of each target is marked at the top of the panel. A downward arrow 
at A re3 t = l-6 or 5/im indicates an upper limit on the stellar or hot dust emission, respectively. 
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NEW RADIO MAPS OF SOURCES IN OUR SAMPLE 
Figure [14] presents the new VLA A-array radio maps. 
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Table 1 

Spitzer HzRG sample and exposure times per instrument. 
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~1 




WN J1123+3141 


11 


23 


55 


74 


+31 


41 


26 


7 


3 


217 


17 


4 


120 


122 


267, 


~1 




PKS 1138-262 


11 


40 


48 


38 


-26 


29 


08 


8 


2 


156 


16 


1 


3000 


122 


9000, -, - 




3C 266 


11 


45 


43 


37 


+49 


46 


08 


2 


1 


275 


17 


1 


120 


305 


267, 


~i 




6C 1232+39 


12 


35 


04 


75 


+39 


25 


38 


9 


3 


220 


18 





120 


122 


267, 


~i 




USS 1243+036 


12 


45 


38 


36 


+03 


23 


20 


7 


3 


570 


19 


2 


120 


122 


267, 


~i 




TN J1338-1942 


13 


38 


26 


09 


-19 


42 


30 


7 


4 


110 


19 


6 


5000 


122 


267, 


~j 




4C 24.28 


13 


48 


14 


87 


+24 


15 


50 


5 


2 


879 






120 


122 


267, 






3C 294.0 


14 


06 


53 


25 


+34 


11 


21 


1 


1 


786 


17 


9 


120 


305 


267, 






USS 1410-001 


14 


13 


15 


10 


-00 


22 


59 


7 


2 


363 






120 


122 


267, 






8C 1435+635 


14 


36 


37 


21 


+63 


19 


14 


4 


4 


250 


19 


4 


120 


122 


134, 


420, 


881 


USS 1558-003 


16 


01 


17 


30 


+00 


28 


46 


2 


2 


527 






120 


122 


267, 






USS 1707+105 


17 


10 


06 


86 


+10 


31 


10 


2 


2 


349 






120 


122 


267, 






LBDS 53W002 


17 


14 


14 


79 


+50 


15 


30 


6 


2 


393 


18 


9 


3300 


122 


134, 


420, 


881 


LBDS 53W091 


17 


22 


32 


93 


+50 


06 


01 


3 


1 


552 


18 


7 


900 


305 


633, 


1311 


, 2643 


3C 356.0 


17 


24 


19 


33 


+50 


57 


36 


2 


1 


079 


16 


8 


120 


305 


134, 


420, 


881 


7C 1751+6809 


17 


50 


50 


03 


+68 


08 


26 


4 


1 


540 


18 


2 


120 


305 


134, 


420, 


881 


7C 1756+6520 


17 


57 


05 


14 


+65 


19 


53 


1 


1 


416 


18 


9 


120 


305 


134, 


420, 


881 


3C 368.0 


18 


05 


06 


37 


+11 


01 


33 


1 


1 


132 


17 


2 


120 


305 


134, 


420, 


881 


7C 1805+6332 


18 


05 


56 


81 


+63 


33 


13 


1 


1 


840 


18 


8 


120 


305 


134, 


420, 


881 


4C 40.36 


18 


10 


55 


70 


+40 


45 


24 





2 


265 


17 


8 


120 


122 


134, 


420, 


881 


TXS J1908+7220 


19 


08 


23 


70 


+72 


20 


11 


8 


3 


530 


16 


5 


120 


122 


134, 


420, 


881 


WIN J1911+634Z 


19 


11 


49 


63 


+63 


42 


09 


6 


3 


590 


19 


9 


120 


122 


134, 


420, 


881 


TN J2007-1316 


20 


07 


53 


26 


-13 


16 


43 


6 


3 


840 


18 


8 


120 


122 


267, 


~ J 




MRC 2025-218 


20 


27 


59 


48 


-21 


40 


56 


9 


2 


630 


18 


5 


120 


122 


267, 






MRC 2048-272 


20 


51 


03 


59 


-27 


03 


02 


5 


2 


060 


18 


3 


120 


122 


267, 






MRC 2104-242 


21 


06 


58 


28 


-24 


05 


09 


1 


2 


491 






120 


122 


267, 






4C 23.56 


21 


07 


14 


80 


+23 


31 


45 





2 


483 


19 


7 


120 


122 


134, 


420, 


881 


MG 2144+1928 


21 


44 


07 


56 


+19 


29 


14 


6 


3 


592 


19 


1 


120 


122 


267, 


~J 




USS 2202+128 


22 


05 


14 


18 


+13 


05 


33 





2 


706 


18 


4 


120 


122 


267, 


~J 




MRC 2224-273 


22 


27 


43 


26 


-27 


05 


01 


7 


1 


679 


18 


5 


120 


305 


267, 






B3 J2330+3927 


23 


30 


24 


82 


+39 


27 


12 


5 


3 


086 


18 


8 


120 


122 


134, 


420, 


881 


4C 28.58 


23 


51 


59 


20 


+29 


10 


29 





2 


891 


18 


7 


120 


122 


267, 






3C 470 


23 


58 


36 


00 


+44 


04 


45 





1 


653 


18 


5 


120 


305 


134, 


420, 


881 



a J2000 coordinates of the identification in the bluest IRAC channel, 
k it-band magnitude within projected 64 kpc radius, taken from S07. 
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Table 2 

Radio data for the SHzRG sample. 





log(L 50 MHz) 


log(L 3 GHz) 




' ' 


^20 


HzRG 


(W Hz" 1 ) 


(W Hz" 1 ) 


Morph a 


(arcsec) 


(%) 


r/" 1 nnoo i /tin 


28.75 


27.75 


rp 
1 


2.3 


2.0 


TV/IT?*" 1 nflQ7 nco 

MHO UUO/-Z00 


07 TO 
Z ( . (A 


OT C\1 

z ( .u ( 


1 


OT R 
Z ( .O 


1 R 
1.0 


OO UUoo+4yo 


27.33 


26.68 


I ~i 
JJ 


3.2 




JVLbvL- U114-zll 


28.66 


28.39 


1 


0.7 


no 

0.3 


UN JUlzl + lozU 


OQ A CI 

zo.4y 


OT A 1 
Z ( .41 


JJ 


n q 




DO UlOz-rOoU 


OT dA 
Z 1 .04 


OR ftfl 

zo.uU 


± 


o o 
Z.Z 


y 


r/^ 1 m /in i one 
OO U14U+OZO 


28.73 


27.89 


JJ 


2.5 




A/fDi^ 1 m co onn 
MHO uioz-zuy 


oq on 
Zo.ZU 


OT TT 
Z ( . / f 


JJ 


o o 
z.Z 


" " " 


MHO UloO-ZDz 


28.46 


27.79 


1 


8.3 


10.0 


ttvt "rnonc i oo/io 
1 IN JUzU0-|-zZ4Z 


OQ A R 


OT A Q 

Z ( .40 


JJ 


O T 
Z. / 




MHO UZll-ZOO 


27.78 


26.22 


c 
o 


<1.5 




TYC HOI 1 1 oo 
IAD UZll-IZZ 


OQ A Q 


OT Q 1 

z ; .ol 


'"T"' 

1 


i t n 
1 ( .U 


Q Q 

. ES 


3C 65 


28.63 


28.06 


1 


1 T C 
1 .0 


0.1 


MHO UZol-ZYo 


28.54 


28.09 


1 


3 . 9 


3 . 5 


MHO UolO-ZOf 


28.95 


28.26 


pv 

JJ 


R T 

0. / 




tv t t~» /~< noo /i o o o 

MHO Uoz4-zzo 


28.49 


27.81 


T 


9.5 


0.2 


MHO 0o50-z79 


28.25 


27.61 


D 


1.2 


' ' ' 


Tv z^ 1 n /I nR O A A 

MHO U4U0-z44 


29.03 


28.11 


1 


10.0 


2.6 


40 bO.07 


29.20 


27.91 


T 


9.0 


2.1 


t> Ty o neon tr , i n 

rKb 0529-549 


29.16 


28.27 


D 


1.2 


<10 


ITS 7"\T Tr\c 1 T i c m o 

WIN J0ol7+5012 


28.02 


26.97 


D 


3.4 




/( /II 1 T 

4o 41.17 


29.18 


28.17 


T 


20.0 


2 . 4 


WIN J0747+3O54 


28.14 


27.02 


S 


2.1 




OOJi, 0ozU-f004z 


28.28 


27.49 


T 


23.4 


0.9 


OO i.ZOi) 


27.82 


27.08 


i -\ 
JJ 


7.6 




Uoo Uozo+iyo 


28.44 


27.47 


T 


12.8 


21.0 


OOH/ uyui+oooi 


OQ 1 O 

zo.iy 


OT A T 
Z I .4 I 


i 


z.y 


n /i 
U.4 


"Do nnno i qa 
JdZ uyuzn~o4 


OQ TQ 


OQ OT 

Zo.Z I 


'"T"' 
1 


1 o 
J.Z 


i q n 


r/^tt 1 nnn c i onct; 

DO-h/ uyUoH-oyoo 


28.17 


27.49 


1 


111 


0.6 


htm Tnno/i oom 
1 IN J(jyz4-ZzUl 


29.51 


27.83 


1 ~\ 
JJ 


1.2 




DO UyOUH-Ooy 


OQ A 1 

Zc5.41 


OT 7fl 

z ( . / y 


1 


1 o 
J.Z 


n q 


TTCC HO/IQ O/IO 

Uoo Uy4o-z4Z 


OQ AO 

Zo.uz 


ot nc 
z ; ,yo 


JJ 


o.y 


" " 


oo zoy 


29.00 


28.19 


1 


11.9 


0.3 


tv^t" 1 imn i ncQ/i 
MO lUiy+UOv54 


OQ CT 
Zf5.0 ( 


OQ 1 Q 

Zo.10 


1 


o o 
Z.Z 


1 T 


MHO 1U1/-ZZU 


OT CiA 

z i .y4 


OQ 1 1 

Zo.ll 


c 
o 


<U.z 




\s/"\t Tn 1 r i rem r 
WIN JlllO+OUlD 


OT QO 
Z ( ,OZ 


OR QT 
ZD.O I 


JJ 


n o 
U.z 




OO zo / 


on i f. 

zy.io 


OQ RO 

Zo.uz 


JJ 


"1 o 

Iz 




\X7"\T T1 1 OQ I Q1/11 

WIN J11Z0+0141 


OQ CI 

Zo.ol 


OT A O 

Z ; .4z 


1 


OC. Q 

zo.o 


tr c 
D.O 


"DT^O 1 1 QQ ORO 

HJ\o lloo-zoz 


29.07 


28.14 


1 


15.8 


3 . 1 


OO ZOO 


OQ C A 

Zo.o4 


ot on 
z ; .oU 


1 ^ 
JJ 


1 R 
J.O 


n i 


R*" 1 1 OQ.O 1 Qfl 

oo izoz+oy 


OQ OQ 

zo.yo 


oq m 
Zo.Ul 


1 


q n 


n t 


TTOC10/1Q I HQR 

Uoolz4o+Uo0 


29.23 


28.25 


1 


6.0 


1 *7 
1. / 


T'TVT T1 QQQ 10/10 

1 IN Jlooo-iy4z 


OQ T1 
Zf5. / 1 


ot nn 
z ( .yu 


1 


C, o 


"1 ("1 

i.y 


A(~^ O/l OQ 

40 Z4.Zo 


on nK 
zy.uo 


OQ OC 

Zo.zo 


1 


O Q 

Z.O 


n t 


o/~i on/i n 

oo zy4.u 


28.96 


28.12 


1 


14.5 


0.3 


ttqq i/iin nn 1 
Uoo 141U-UU1 


28.41 


27.69 


1 


24.0 


6.7 


oO 14OO+D00 


on a n 

zy.4u 


OQ C C 

zo.oo 


1 


4.0 


. 


TTQQ 1 CCO AAO 

Uoo looo-UUo 


28.82 


28.00 


1 


9.2 


1.6 


TTQQ 1 '7CV7 1 1 AC 


28.63 


27.78 


JJ 


22.5 




LJoJJo OoWUUz 


27.78 


27.02 


o 
o 


<1.5 




LUDb ooWUyl 


27.04 


26.29 


D 


4.3 


' ' ' 


oa oce n 
OO ooO.U 


28.35 


27.65 


i "T"l 
1 


76.2 


0.1 


T/^ 1 1 TCI 1 COfifl 

/o 1 ( oi+ooUy 


27.46 


27.01 


1 ~\ 
JJ 


2.0 




/O 1 ( OD+DOZU 


27.40 


27.00 


JJ 


O T 
Z. / 


' ' ' 


OO ODO.U 


28.52 


27.63 


1 


8.8 


0.2 


7o 18U5+boOZ 


27.78 


27.12 


T 


15.6 


0.6 


/i (~* /in 

40 4U.OD 


OQ 7fl 

zc5. / y 


ot n/i 
z ; .y4 


JJ 


i n 




rpv q t 1 nno i to on 
lAo JlyUo+/zzU 


29.12 


28.15 


1 


15.4 


4.8 


WN J1911+6342 


28.14 


27.03 


s 


1.8 




TN J2007-1316 


29.13 


27.79 


T 


3.5 


15 


MRC 2025-218 


28.74 


27.96 


T 


5.1 


0.7 


MRC 2048-272 


28.72 


27.85 


D 


6.7 


<0.2 


MRC 2104-242 


28.84 


27.88 


T 


23.7 


0.7 


4C 23.56 


28.93 


28.26 


T 


53.0 


14.6 


MG 21444-1928 


29.08 


28.27 


D 


8.9 




USS 2202+128 


28.54 


27.75 


T 


4.2 


0.9 


MRC 2224-273 


27.52 


27.39 


S 


<0.2 




B3 J2330+3927 


28.33 


27.59 


T 


3.0 


50 


4C 28.58 


28.91 


27.89 


T 


14.5 


2.0 


3C 470 


28.79 


28.24 


T 


24.6 


0.4 


Radio source morphol 


ogy: S = Singlc compoi 


cnt, D=Double co 


mponent, T= 


= Three or m 




Largest angular size c 


f the radio source; see 


Carilli et al. \ 1997). 






Core fraction at reatfi 


ame 20 GHz; seeldaril 


li et al.| I1997IV 
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The Spitzer High Redshift Radio Galaxy Survey 



Table 3 

Spitzer photometry for HzRGs (fijy). 



HzRG 



/3.6/im 



/4.5/im 



/5.8/im 



/s.o/im 



fwfira /24/^m 



/l60/im 



6C 0032+412 
MR C0037-258 
6C*0058+495 
MRC 0114-211 
TN J0121+1320 
6C* 0132+330 
6C 0140+326 
MRC 0152-209 
MRC 0156-252 
TN J0205+2242 
MRC 0211-256 
TXS 0211-122 

3C 65 
MRC 0251-273 
MRC 0316-257 
MRC 0324-228 
MRC 0350-279 
MRC 0406-244 

4C60.07 
PKS 0529-549 
WN J0617+5012 

4C41.17 
WN J0747+3654 
6CE 0820+3642 

5C7.269 
USS 0828+193 
6CE 0901+3551 

B2 0902+34 
6CE 0905+3955 
TN J0924-2201 
6C 0930+389 
USS 0943-242 

3C 239 
MG 1019+0534 
MRC 1017-220 
WN J1115+5016 

3C 257 
WN J1123+3141 
PKS 1138-262 

3C 266 
6C 1232+39 
USS 1243+036 
TN J1338-1942 
4C 24.28 
3C 294.0 
USS 1410-001 
8C 1435+635 
USS 1558-003 
USS 1707+105 
LBDS 53W002 
LBDS 53W091 

3C 356a 
7C 1751+6809 
7C 1756+6520 

3C 368.0 
7C 1805+6332 
4C 40.36 
TXS J1908+7220 
WN J1911+6342 
TN J2007-1316 
MRC 2025-218 
MRC 2048-272 
MRC 2104-242 

4C 23.56 
MG 2144+1928 
USS 2202+128 
MRC 2224-273 
B3 J2330+3927 
4C 28.58 
3C 470 



15.4+ 2.1 
221.0+22.0 
82.0+ 8.3 
87.3+ 8.9 
10.5+ 1.8 
31.6+ 3.6 

< 623.0 
108.0+11.0 
291.0+29.0 

7.5+ 1.7 

166.0+17.0 
30.2+0.7 

161.0+32.0 
11.5+ 1.6 
19.3+ 2.1 
39.4+ 4.2 
23.6+ 2.6 
40.4+ 4.3 
20.6+ 3.1 
46.6+ 4.9 
3.6+ 1.0 
23.4+ 2.4 
19.1+ 2.4 
79.2+ 8.1 
41.0+ 4.5 
61.7+ 6.9 
37.2+ 4.1 
6.4+ 0.8 
51.8+ 5.4 
11.3+ 1.8 
30.7+ 3.4 
21.5+ 2.6 
96.4+ 9.8 
25.6+ 2.9 

119.0+12.0 
7.8+ 2.1 
85.0+ 8.7 
48.2+ 5.0 

318.0+32.0 
67.9+ 7.0 
33.3+ 3.6 
22.0+ 2.5 
17.8+ 1.9 
16.5+ 2.2 

< 93.0 
50.6+ 5.3 
18.6+ 2.1 
78.8+ 8.1 
22.1+ 2.7 
32.0+ 3.3 
43.3+ 6.3 
108.0+11.0 
46.6+ 4.9 
39.6+ 4.2 

126.0+13.0 
28.4+ 3.6 
36.5+ 3.9 
200.0+20.0 
9.5+ 1.8 
56.2+ 5.9 
68.4+ 7.1 
59.5+ 6.2 
28.1+ 3.3 
61.1+ 6.4 
22.1+ 2.7 
60.4+ 6.3 
61.6+ 6.4 
99.6+10.1 
31.6+ 3.5 
49.5+10.4 



33.8+ 3.7 
248.0+25.0 

86.7+ 8.8 
117.0+12.0 

14.4+ 2.1 

41.3+ 4.5 

< 450.0 
165.0+17.0 
405.0+41.0 

4.7+ 2.3 
197.0+20.0 
47.4+1.1 

12.7+ 1.8 
20.1+ 2.1 
39.7+ 4.3 
40.6+ 4.2 
43.3+ 4.6 
27.0+ 2.0 
52.9+ 5.5 
5.5+ 1.2 
27.5+ 2.8 
25.3+ 3.0 
81.9+ 8.4 
49.5+ 5.3 

133.0+13.0 
46.5+ 5.0 
9.9+ 1.1 
60.1+ 6.2 
11.4+ 1.9 
32.2+ 3.6 
28.4+ 3.2 

111.0+11.0 
19.5+ 3.8 

179.0+18.0 
9.5+ 3.0 

111.0+11.0 
74.4+ 7.6 

497.0+50.0 
73.1+ 7.5 
41.8+ 4.4 
21.5+ 3.0 
10.7+ 1.2 
27.3+ 3.1 

< 103.0 
79.0+ 8.1 
14.7+ 2.9 

101.0+10.3 
30.1+ 3.4 
44.0+ 4.5 
51.4+ 5.3 

110.0+11.0 
50.8+ 5.3 
46.9+ 5.0 

112.0+11.0 
42.1+ 4.4 
41.3+ 4.3 

229.0+23.0 
8.7+ 1.8 
54.4+ 5.8 
77.1+ 8.0 
72.6+ 7.5 
29.7+ 3.5 
86.2+ 8.8 
18.3+ 2.4 
95.8+ 9.8 
86.1+ 8.8 

143.0+14.0 
36.0+ 3.9 
75.2+11.8 



57.6+ 8.9 
286.0+29.0 

93.2+ 9.5 
157.0+16.0 

< 32.5 
51.9+12.4 

< 34.9 
215.0+22.0 
717.0+72.0 

< 29.3 
222.0+22.0 

126.5+6.3 
320.0+66.0 

< 29.1 
19.5+ 2.4 
61.1+ 8.6 
82.2+27.6 

< 51.6 

< 26.0 
62.7+ 8.6 

< 48.3 
35.6+ 3.7 

< 29.9 
82.0+ 8.4 
57.8+ 6.1 

201.0+21.0 
52.8+10.7 
11.0+ 2.5 
96.8+ 9.8 

< 30.6 
37.3+ 9.2 

< 30.9 
130.0+12.0 

< 35.4 
273.0+27.0 

< 54.8 
194.0+19.0 

92.7+ 9.4 
887.0+89.0 
45.1+ 4.7 
52.0+ 9.2 

< 48.0 
14.3+ 2.8 
43.3+ 8.3 
68.0+16.8 
166.0+17.0 

< 45.5 
105.0+10.5 

22.6+ 9.1 
49.9+ 5.2 
23.9+ 6.1 
122.0+14.0 

< 40.9 
34.7+ 7.7 

112.0+11.0 
51.4+ 5.4 
45.4+12.9 

241.0+25.0 

< 26.4 
41.1+10.7 
86.8+10.9 
78.3+10.4 
32.8+10.0 
126.9+12.8 

< 25.5 
120.1+12.1 
98.4+10.0 
160.0+16.0 

< 41.7 
70.9+10.4 



98.3+ 10.0 
518.0+ 52.0 
309.0+ 31.0 
398.0+ 40.0 

< 40.3 
109.0+ 11.0 

<36.3 
415.0+ 42.0 
1125.0+113.0 

< 37.4 
278.0+ 28.0 

462.2+7.9 

< 28.8 
38.1+ 4.1 
89.9+ 9.9 
79.3+ 26.5 
63.5+ 14.5 

< 39.0 
72.2+ 9.1 

< 54.5 
36.5+ 3.5 
44.9+ 11.5 
68.0+ 7.0 

< 40.1 
687.0+ 74.0 
69.8+ 12.4 
41.3+ 2.3 

146.0+ 14.0 

< 33.0 

< 32.6 
25.8+ 11.7 
142.0+ 14.0 

< 42.9 
360.0+ 36.0 

< 61.1 
322.0+ 33.0 
182.6+ 18.4 

1500.0+150.0 
102.6+ 10.4 
75.2+ 10.8 

< 62.9 
9.9+ 3.4 

102.0+ 10.0 
66.6+ 20.6 
240.0+ 24.0 

< 50.3 
233.0+ 23.4 

< 33.3 
103.0+ 11.0 

26.6+ 6.4 
434.0+ 47.0 
36.5+ 16.0 
42.6+ 8.6 
210.0+ 21.0 
95.6+ 17.1 
26.3+ 10.3 
480.0+ 48.0 

< 26.0 
121.3+ 12.3 
126.8+ 12.9 
38.4+ 13.1 

< 36.3 
423.7+ 42.5 

< 30.3 
106.6+ 10.8 
203.9+ 20.5 
474.0+ 47.0 

40.9+ 4.4 
266.0+ 30.0 



412+ 83 
877+100 
1280+100 
1690+100 
<204 
<280 
<160 
1580+100 
1980+100 

<211 
378+ 86 
1590+220 
1220+330 
102+ 33 

<156 
530+ 54 
208+ 40 
637+ 86 

a 

248+ 89 
<203 
<181 
<198 
225+ 70 

<260 
1910+130 
410+ 70 
450+100 
1440+120 
<168 
<180 
170+ 48 
848+ 90 
220+ 70 
740+120 

<160 
1280+150 
1160+110 
3020+100 
660+120 
360+ 80 
<260 
<226 
573+ 96 
430+ 99 
660+120 

<86 
993+100 

<123 
593+114 

<166 
2280+220 
147+ 57 
220 a ±114 
1620+180 
310+ 80 
54+ 28 
1410+170 

<137 
378+113 
200+ 62 
<220 
<217 
2400+ 90 

<244 
417+ 88 
625+117 
1170+ 90 
430+ 85 
1540+180 



467+ 42 
1740+ 39 
1490+ 50 
2090+ 40 

<131 
109+ 25 
<2190 
3320+133 

<116 
710+ 35 
2750+ 40 
2072+ 30 
476+ 33 

<102 
1880+ 35 
306+ 44 
1540+ 40 

a 

966+ 40 
<122 

370+ 40 
<117 

1290+ 40 

2880+ 40 

1340+ 40 
323+ 12 

3700+ 40 

<125 
530+ 43 
518+ 40 

1890+ 60 
415+ 43 

1140+ 30 
<100 

1770+130 

1670+ 37 

3890+ 20 
985+ 22 
459+ 33 
511+ 43 
<99 

1140+ 40 
471+ 31 

1070+ 40 
<119 

1270+ 37 
264+ 32 
648+ 40 
<45 

4170+ 40 
342+ 34 
95 a ±40 

3350+ 50 
673+ 35 
536+ 37 

1910+ 49 

<102 
385+ 40 
216+ 43 

a 

709+ 48 
4630+ 40 

529+ 33 
1341+ 35 
1060+ 40 
2320+ 40 

866+155 
2750+ 40 



< 3950 
18900+1962 



2580 



< 1040 
24700+2306 

< 3750 

< 4110 

< 3210 



3330 
3390 



4130 



< 4300 

< 1610 

< 4400 

< 3600 

< 5900 
28800+2710 

< 4330 

< 4750 
16200+1905 

< 3710 



< 96000 

< 94600 

< 20700 



< 88800 

< 47700 

< 64600 

< 74100 

< 26300 



< 18500 

< 52900 

< 50900 



< 33900 



< 65100 

< 23900 

< 70200 

< 51600 

< 151000 

< 39000 

< 111000 

< 64500 

< 63300 

< 57100 



30300+2958 < 70500 



< 4670 

< 5570 



< 64300 

< 102000 



rby object. 
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Table 4 

Results of SED fitting. 



Source 


log(Li.6 fim/L©) 


log(i3Atm/L ) 


log(Ls /Lim/L©) 


fa 


log(^H,*/L©) 


log(M*/M ) 


6C 0032+412 


12.07 


12.28 


12.29 


0.21 


<11.37 


<11.18 


MRC 0037-258 


11.49 


11.50 


11.53 


0.95 


11.44 


11.56 


6C* 0058+495 


11.19 


11.20 


11.53 


0.99 


11.16 


11.26 


MRC 0114-211 


11.40 


11.70 


11.96 


0.96 


11.36 


11.39 


TN J0121+1320 


11.55 


11.36 


11.73 


0.46 


11.19 


11.02 


6C* 0132+330 


11.14 


11.33 


11.23 


0.84 


<11.04 


<11.03 


6C 0140+326 


11.84 


12.61 


12.89 


0.88 


<11.76 


<11.42 


MRC 0152-209 


11.88 


12.15 


12.26 


0.91 


11.80 


11.76 


MRC 0156-252 


12.40 


12.57 


12.43 


0.54 


12.11 


12.05 


TN J0205+2242 


11.03 


11.42 


11.79 


0.97 


10.99 


10.82 


MRC 0211-256 


11.58 


11.49 


11.24 


0.83 


<11.47 


<11.54 


TXS 0211-122 


11.80 


12.43 


12.40 


0.31 


<11.26 


<11.16 


3C 65 


11.43 


11.66 


11.73 


0.92 


<11.37 


<11.47 


MRC 0251-273 


11.33 


11.57 


11.88 


0.66 


11.12 


10.96 


MRC 0316-257 


11.49 


11.77 


11.70 


0.79 


11.37 


11.20 


MRC 0324-228 


11.33 


11.44 


11.71 


0.96 


11.29 


11.25 


MRC 0350-279 


11.24 


11.44 


11.30 


0.69 


<11.04 


<11.00 


MRC 0406-244 


11.53 


11.65 


12.10 


0.99 


11.49 


11.38 


4C 60.07 


11.73 


12.51 


12.91 


0.88 


11.65 


11.44 


PKS 0529-549 


11.66 


11.75 


11.95 


0.94 


11.60 


11.46 


WN J0617+5012 


11.18 


11.38 


11.49 


0.37 


10.72 


10.55 


4C 41.17 


11.66 


12.00 


12.21 


0.93 


11.60 


11.39 


WN J0747+3654 


11.55 


11.43 


11.38 


0.63 


11.33 


11.15 


6CE 0820+3642 


11.57 


11.19 


11.27 


0.95 


11.52 


11.48 


5C 7.269 


11.51 


11.06 


11.89 


0.95 


11.46 


11.38 


USS 0828+193 


12.18 


12.75 


12.68 


0.40 


<11.75 


<11.60 


6CE 0901+3551 


11.34 


11.37 


11.61 


0.97 


11.30 


11.25 


B2 0902+34 


11.52 


12.15 


12.06 


0.30 


<10.97 


<10.81 


6CE 0905+3955 


11.47 


11.73 


12.16 


0.98 


11.43 


11.39 


TN J0924-2201 


12.19 


12.40 


12.18 


0.21 


11.50 


11.10 


6C 0930+389 


11.32 


11.44 


11.80 


0.99 


11.28 


11.17 


USS 0943-242 


11.43 


11.49 


11.93 


0.99 


11.40 


11.22 


3C 239 


11.66 


11.57 


11.90 


0.99 


11.63 


11.60 


MG 1019+0534 


11.36 


11.53 


11.89 


0.99 


11.32 


11.15 


MRC 1017-220 


11.84 


11.97 


11.85 


0.81 


<11.72 


<11.70 


WN J1115+5016 


10.95 


11.48 


11.55 


0.82 


10.84 


10.70 


3C 257 


12.03 


12.42 


12.44 


0.81 


<11.91 


<11.78 


WN J1123+3141 


12.13 


12.67 


12.68 


0.60 


<11.89 


<11.72 


PKS 1138-262 


12.52 


12.84 


12.73 


0.68 


<12.33 


<12.26 


3C 266 


11.13 


10.90 


11.10 


0.99 


11.10 


11.18 


6C 1232+39 


11.80 


12.18 


12.15 


0.73 


<11.64 


<11.47 


USS 1243+036 


11.62 


12.17 


12.28 


0.71 


<11.45 


<11.27 


TN J1338-1942 


11.37 


11.55 


11.92 


0.99 


11.34 


11.04 


4C 24.28 


11.58 


12.18 


12.27 


0.53 


<11.28 


<11.11 


3C 294.0 


11.41 


11.29 


11.56 


0.99 


11.38 


11.36 


USS 1410-001 


11.92 


12.02 


11.86 


0.42 


<11.52 


<11.41 


8C 1435+635 


11.80 


12.04 


12.08 


0.56 


11.53 


11.21 


USS 1558-003 


11.95 


12.05 


12.05 


0.82 


<11.83 


<11.70 


USS 1707+105 


11.29 


10.86 


10.96 


0.89 


11.21 


11.11 


LBDS 53w002 


11.55 


11.72 


11.68 


0.71 


<11.38 


<11.27 


LBDS 53w091 


11.20 


10.59 


10.20 


1.00 


11.18 


11.19 


3C 356.0 


11.21 


11.23 


11.62 


0.99 


11.18 


11.29 


7C 1751+6809 


11.20 


10.78 


10.91 


1.00 


11.17 


11.19 


7C 1756+6520 


11.15 


10.92 


11.32 


1.00 


11.12 


11.15 


3C 368.0 


11.35 


11.10 


11.41 


1.00 


11.32 


11.43 


7C 1805+6332 


11.21 


11.39 


11.40 


0.84 


<11.11 


<11.07 


4C 40.36 


11.43 


10.93 


11.14 


0.95 


11.38 


11.29 


TX J1908+7220 


12.70 


12.91 


12.82 


0.59 


<12.44 


<12.27 


WN J1911+6342 


11.38 


11.94 


11.87 


0.57 


11.11 


10.93 


TN J2007-1316 


12.19 


12.42 


12.24 


0.56 


<11.92 


<11.69 


MRC 2025-218 


11.87 


11.86 


11.69 


0.87 


<11.78 


<11.62 


MRC 2048-272 


11.56 


11.30 


11.68 


1.00 


11.53 


11.47 


MRC 2104-242 


11.35 


11.55 


11.94 


0.98 


11.32 


11.19 


4C 23.56 


11.89 


12.60 


12.74 


0.71 


<11.71 


<11.59 


MG 2144+1928 


11.49 


12.22 


12.32 


0.70 


<11.31 


<11.13 


USS 2202+128 


11.97 


11.61 


11.89 


0.70 


11.78 


11.62 


MRC 2224-273 


11.47 


11.63 


11.78 


0.95 


11.42 


11.41 


B3 J2330+3927 


12.43 


12.71 


12.66 


0.51 


<12.11 


<11.94 


4C 28.58 


11.57 


11.84 


12.24 


0.98 


11.54 


11.36 


3C 470 


11.35 


11.72 


12.02 


0.95 


11.30 


11.30 



a Stellar fraction at 1.6 fim in our SED model (see §3.2). 
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Table 5 

Binary AGNs 



Source 


RA 
(J2000) 


Dec. 
(J2000) 


fa.6 fim 

(**Jy) 


74.5 

(**Jy) 


/5.8/im 

(/xJy) 


fs.O )im 


/16 fim 

(^Jy) 


724 jim 

(^Jy) 


4C 60.07B 
3C 3566 
CI 1756.7 
MRC 2048- 272B 


05 12 54.76 
17 24 19.02 
17 57 04.98 
20 51 03.30 


+60 30 48.9 
+50 57 40.3 
+65 19 51.0 
-27 03 03.2 


<3.3 
72.6+7.3 
23.0+2.3 
24.6+2.5 


7.8+1.0 
71.1+7.2 
33.7+3.4 
32.9+3.4 


23.5+3.3 
69.4+7.5 
28.1+3.8 
35.2+4.7 


43.3+ 5.1 
231.9+23.4 
37.6+ 4.7 
34.3+ 5.0 


555+ 83 
2280+220 
220+150 
<220 


1390+50 
4170+40 
370+40 
491+42 



